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Nonlinear interaction of spin and charge currents in graphene 
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We describe a nonlinear interaction between charge currents and spin currents which arises from 
the energy dependence of the conductivity. This allows nonmagnetic contacts to be used for measur- 
ing and controlling spin signals. We choose graphene as a model system to study these effects and 
predict its magnitudes in nonlocal spin valve devices. The ambipolar behavior of graphene is used 
to demonstrate amplification of spin accumulation in p-n junctions by applying a charge current 
through nonmagnetic contacts. 

PACS numbers: 72.25.Hg, 72.80. Vp, 75.76.+j, 85.75.-d 



Spin-polarized transport (spintronics) [1] in graphene, 
a onc-atom-thick layer of carbon [2], is of both funda- 
mental and technological interest due to its long spin 
relaxation length [3] and large spin signals [4]. In this 
work we focus on understanding graphene spintronics as 
it is experimentally addressed by an all-electrical scheme 
involving the use of a nonlocal device geometry [3]. The 
conductivity in graphene has been considered only at the 
Fermi level, which leads to equal values for both spin 
channels. We point out that very recent work [5] uses the 
energy dependent conductivity of graphene, in conjunc- 
tion with Zccman splitting by applied magnetic fields, to 
explain a giant spin-Hall effect in the linear regime. 

In this contribution we highlight nonlinear effects in 
the absence of external magnetic field that gives rise to in- 
teractions between electron spin and charge in graphene. 
We argue that the nonlinear interaction between spin and 
charge can be consistently described by using the equa- 
tions for spin diffusion in metals [6, 7], while considering 
the energy dependent conductivity a{t) of graphene and 
the large spin accumulation achievable by spin injection. 
This gives rise to phenomena observable in the nonlocal 
geometry, which include detection of spin accumulation 
in graphene via nonmagnetic contacts, its manipulation 
using charge currents and amplification in bipolar de- 
vices. 

Previous experimental work has shown the manipula- 
tion of spin accumulation in graphene by applying high 
electric fields [8] . Such a manipulation has been later in- 
terpreted by considering the effect of low-resistance con- 
tacts [9] within the drift-diffusion formalism for spin ac- 
cumulation derived for semiconductors [10]. In the fol- 
lowing, we consider highly resistive noninvasive contacts 
[11, 12], which can be treated as ideal (spin) current in- 
jectors and (spin) voltage detectors [4, 13]. This allows 
us to only focus on the physics of spin transport within 
graphene. 

We start with the well established model for spin trans- 
port in metals [6, 7] with the electrochemical potential for 
each spin channel expressed as = /Xavg ± A//, where 
the index ± refers to electron spin ±1/2 . Here, A// is 
the term related to the spin accumulation and /.tavg the 



average potential. This results in spin-diffusion equations 
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with e the electron charge, A the spin relaxation length, 
a± and J± the conductivities and the currents for each 
spin channel, and E{x) — {l/e){d^a,vg/dx). As shown in 
Eq. (2), the gradient in A^ drives the current for each 
spin channel in opposite directions, whereas the electric 
field E drives them in the same direction. The spin- 
dependent conductivities are described by the polariza- 
tion /3 such that a± = [2p(l ± /3)]~^. The general solu- 
tions for fi± , J± , and E in a, homogeneous medium were 
presented in Appendix C of Ref. 7. To find numerical 
solutions for a certain device configuration we divide the 
graphene into discrete regions and use the solutions for 
each region, while keeping continuities of /x± and J± . 

In metals it is possible to achieve a spin accumula- 
tion A/i « 10 /ieV [14] whereas in graphene it can be 
considerably larger, as explained below. Electrical spin 
injection from a tunnel contact with polarization P using 
a charge current Ii results in injection of a spin current 
PIi- Ignoring any interaction between spin and charge, 
Afj, decays away from the injector as given by [3] 
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with p = 1/(7 the graphene square resistance, w its width, 
and Rs the nonlocal spin resistance as measured by a sec- 
ond magnetic contact with same P. For typical values of 
p 1 kri, w = 1 /xm, A = 2 fim, Ii = 10 /lA and 
P = 20 % the resulting A/i profile reaches ~ 1 mcV, 
as depicted in Fig. 1(a). At such large splitting in 
the electrochemical potential the energy dependence of 
the graphene conductivity a{e) starts to be noticeable. 
Therefore we introduce a splitting dependent /3 given by 
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FIG. 1. Generation of a nonlocal charge voltage Vni ~ 
— /iavg/e by spin injection in graphene. (a) Profiles of /iavg 
and A/i created by a spin current. The linear ohmic drop due 
to Ji on the right side of the circuit has been subtracted from 
/iavg for the sake of clarity, (b) Charge carrier density depen- 
dence of the nonlocal potential 1 fim away from the injector. 
The inset shows Vni measured by a nonmagnetic contact ver- 
sus Ii. Dashed curves correspond to injection of pure spin 
current PIi without a charge current, whereas solid curves 
are for considering the effect of a charge current h in the 
right side of the circuit. 



with a conductivity spin polarization proportional to A/i 
and a coefficient a. The effects of temperature and dis- 
order, described later, can be taken into account by con- 
sidering the experimental electrical conductivity cr(e,T). 

Now, we discuss the nature of the coefficient a. 
The conductivity of graphene away from the neutrality 
(Dirac) point can be described by ct = vne, with v the 
carrier mobility and n the carrier density. Due to the 
linear density of states in graphene [2] the carrier den- 
sity depends on energy as n(e) = e^/nfi^Vp, where ft is 
the reduced Planck constant and vp ~ 10^ m/s is the 
Fermi velocity. For a constant mobility, the latter leads 
to a = 2/e cx l/y/n, so a diverges when n tends to zero. 
Therefore the maximum value of a is given by the mech- 
anisms limiting how close we can get to n = 0. The 



charge carrier density in graphene field-effect transistors 
is electrostatically tunable by applying a voltage Vg to a 
gate such that Ug = {Cg/e){Vg — Vd) with Cg the gate 
capacitance per unit area and Vd the voltage at which the 
neutrality point is located. A background carrier density 
Tii is induced by the presence of electron-hole puddles 
(upd) due to disorder [15] and thermal generation of car- 
riers (nth) [16]. Up to date, all spintronic devices have 
been fabricated with graphene supported on a substrate, 
mostly Si02. At room temperature, for a typical value 
of m = {Up^ + 4ri,j^)^/^ = 4 X 10^-^ cm~^, we obtain a 
maximum value of a w (60 mcV)^^ which together with 
A/i « 1 meV yields /3 « 1%. 

We use the modeling above to study spin transport in 
graphene field-effect transistors. We can describe most 
experimental a vs Vg (Dirac) curves for graphene on 
Si02 by taking a constant value oi v = 0.4 m^/Vs and 
a carrier density n = {rig + nfy/^ [16]. For simplicity, 
we keep the polarization of magnetic contacts fixed at 
P = 20 % and a carrier density independent spin relax- 
ation length of A = 2 /tm. Unless stated otherwise, we 
use Ug ~ 2 X 10^^ cm~^ [a w (83 mcV)~^], well into the 
metallic regime. First, we consider in Fig. 1(a) a non- 
local measurement where a spin current PIi is injected 
into graphene. The presence of a spin current in graphene 
with an inhomogeneous conductivity polarization /3 cre- 
ates a nonlocal charge voltage V^i = — /iavg/e- Interest- 
ingly, such a potential is detectable with both magnetic 
and nonmagnetic contacts. For the simple case of pure 
spin current injection (ignoring the effect of the charge 
current Ii on the right side of the circuit) we obtain 
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which indicates that this is a second order effect in A/t. 
The latter also means that Vni must decay with a char- 
acteristic length of A/2. The nonlinear nature of Vni is 
explicit in the inset of Fig. 1(b). This effect opens the 
unique possibility of measuring spin signals in graphene 
without using magnetic detectors. 

Next, we consider the effect of changing the graphene 
carrier density Ug on V^i- We calculate a from 
Eq. (4) using the simulated electrical conductivity 
a{e, T), which includes the effect of rn, such that da/de = 
{da / dng){dng / de) . The coefficient a has opposite po- 
larity for hole and electron regimes. Besides, a finite 
Tii introduces the presence of both electrons and holes 
near the Dirac point, which removes the divergence of a 
and makes it zero. This behavior is in analogy to the 
case of the thermoelectric Seebeck coefficient S which 
has the same dependence on tj(e) as a and follows an ap- 
proximate Mott formula based on the experimental Dirac 
curve [17]. From Eqs. (3), (4) and (5) the dependence of 
Vn\ on energy is given by a~^da/de. 

The nonlocal charge voltage goes as Vn\ oc if for injec- 
tion of pure spin currents. If we also consider the charge 



current /i on the right side of the circuit [Fig. 1(a)], 
an asymmetry in the Vni vs Ii curve is visible [inset in 
Fig. 1(b)]. This is a higher (> 3"^) order effect on Vni 
given by the interaction of Ii with the spin accumulation 
on the right side of the circuit, which creates a higher 
(> 2"^) order effect on the nonlocal spin accumulation 
A/x. We make explicit such an effect on A// in Fig. 2(a). 
The resulting nonlocal spin resistance detected in a spin 
valve with a second magnetic contact, Rs = PA/i/e/i, 
now varies with the injection current, as observed in the 
nonlinear behavior shown in the inset of Fig. 2(a). 

An interesting result is observed if we consider a second 
charge current I2 via a nonmagnetic contact, as depicted 
in Fig. 2(b). A spin accumulation in graphene creates a 
conductivity spin polarization /3, which in the presence of 
a charge current I2 gives rise to a spin current /3/2. The 
nonmagnetic contact hence seems to inject a spin current 
similarly to the case of a magnetic contact. Depending on 
the polarities of I2 and a, spin accumulation or depletion 
is observed. The nonmagnetic contact offers an extra 
handle by which we can amplify the spin accumulation. 
As shown in Fig. 2(b), A^ under the nonmagnetic contact 
can be even larger than that under the magnetic contact. 

In graphene field-effect devices we can individually ad- 
dress specific regions via local electrostatic gates. We 
use this capability to study ambipolar spin transport in 
graphene. We choose a highly symmetric case where the 
physics can be easily understood and derive a simple an- 
alytical description which accurately describe the simula- 
tions. The latter is possible because in graphene we can 
ignore the effects of the charge depletion region present 
in nondegenerate p-n junctions [1, 18]. 

We consider a graphene channel with the left half set 
in the hole regime and the right half set in the electron 
regime, as depicted in Fig. 3. Initially, with I2 = 0, the 
magnetic contact located at the junction creates a spin 
accumulation A^o given by Eq. (3) (assuming pure spin 
current injection). We define A = Xpe/2w, so that under 
the magnetic contact we have an initial A/xq = APIi. 
If we now apply a charge current I2 via the nonmag- 
netic contacts, the sign change of the parameter a at 
the junction creates a source of spin current equal to 
2/3/2. Such a discontinuity induces a spin accumulation 
A/Xind- We remark that the graphene spin polarization 
/3 is given by the total spin accumulation at the junction 
A/itot = A/iQ + A/iind- Therefore at the magnetic contact 



A/xtot = APh 
^ = A2al2 



A2al2A^tot + S. 
9A^tot(a;) 
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dx , (7) 



with ^ a compensation term in A/Xind corresponding to 
the presence of I2 in the p and n regions with inhomo- 
geneous spin polarization /?. For small charge currents 
{I2 ^ 1/aA) we have A/xtot(a^) ~ A/xo(x) and the inte- 
gral in Eq. (7) evaluates to — A/itot/2. Introducing this 
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FIG. 2. Effect of charge current on spin accumulation Afi. (a) 
Profile of Afj. due to its interaction with Ii on the right side 
of the circuit. Inset: nonlinearity of the nonlocal A/x, 1 fim 
away from the injector, (b) Redistribution of A/x caused by 
a current I2 via a nonmagnetic contact. 



result into Eq. (6) leads to A/xtot at the junction 
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equivalent to an amplifier circuit with positive feedback 
controlled by a/2. Eq. (8) gives accurate results for low 
values of h- The divergence at Aah = 1 is a result 
of our approximation for ^. In reality, the distribution 
of spin accumulation will (de)focus at the junction with 
changing I2 [18] which yields different compensation ^. 

To account for large values of I2 we generalize Eqs. (6) 
and (7) for the case of spin accumulation A/xtot(x) at any 
location within the graphene. It follows that the general 
solution of A/itot satisfies 
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where the second term arises due to the spin-dependent 
conductivity. We describe electronic transport in energy 
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FIG. 3. Change in spin accumulation profile in a graphene 
p-n junction due to a charge current I2 (assuming pure spin 
current injection P7i). Inset: amplification of at the junc- 
tion (under the magnetic contact) due to 



space via a. For a one-dimensional Drude model the 
mathematical formulation is similar to that of drift [10]. 

Eq. (9) has solutions of the form exp(=Fa::/i±) with 
L± = \/{±Aal2 + ^^{AahY + !)• Using these solu- 
tions together with Eqs. (6) and (7) we find that, for the 
case of nonmagnetic contacts far away from the spin in- 
jector, the spin accumulation has the form A/itot(a;) = 
(AP/i/A)L_ exp (— The analytical solution de- 
scribes the (de)focusing of the A/i profile with I2 shown 
in Fig. 3. At /2 > the distribution of A/i focuses at the 
junction. The opposite occurs for I2 < 0. In the limit 
/2 3> the peak in the spin accumulation has a value 
of A/x = Pa/1/2 and the distribution tends towards 
a Dirac delta function with constant area APIiX. The 
small peaks at a; = ±9 /im arc due to amplification at 
the nonmagnetic contacts. The dependence of A/i at the 
junction versus I2 is shown in the inset of Fig. 3. Fi- 
nally, note that using the solutions to Eq. (9) we can 
also describe the nonlinear spin resistance caused by Ii 
in Fig. 2(a), as A/i = (PpAeffe/i)(2ui)~-^ exp(— |x|/A) for 
a; < 0, where Aoff = 2(1/A + l/L-)'^. For small h the 
latter leads to the addition of a second order term to 
Eq. (3) of the form Kie/P, with Ki defined in Eq. (5). 

In conclusion, we have described the interaction be- 
tween spin and charge transport in graphene by treating 
it as a ferromagnet with a conductivity spin polariza- 
tion P induced by the presence of a spin accumulation 
A/t. This leads to phenomena experimentally accessible 
via nonlocal measurements, including detection and ma- 
nipulation of spin signals with nonmagnetic contacts, its 
dependence on carrier density and amplification effects in 



ambipolar devices. Since the nonlinear interaction arises 
solely due to the energy dependence of the conductivity, 
the ideas described in this work are also applicable to 
other materials used for spin transport, such as Si and 
GaAs. The generality of this interaction is analogous 
to the interaction between heat and charge described by 
thermoelectricity. 
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